Traffic properties for stochastic routings on scale-free networks 
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For realistic scale-free networks, we investigate the traffic properties of stochastic routing inspired 
by a zero-range process known in statistical physics. By parameters a and 5, this model controls 
degree-dependent hopping of packets and forwarding of packets with higher performance at more 
busy nodes. Through a theoretical analysis and numerical simulations, we derive the condition for 
the concentration of packets at a few hubs. In particular, we show that the optimal a and & are 
involved in the trade-off between a detour path for a < and long wait at hubs for a > 0; In the 
low-performance regime at a small 5, the wandering path for q < better reduces the mean travel 
time of a packet with high reachability. Although, in the high-performance regime at a large 5, 
the difference between a > and a < is small, neither the wandering long path with short wait 
trapped at nodes (a = — 1), nor the short hopping path with long wait trapped at hubs (a = 1) is 
advisable. A uniformly random walk (a = 0) yields slightly better performance. We also discuss the 
congestion phenomena in a more complicated situation with packet generation at each time step. 

PACS numbers: 89. 75. He, 02.50.Ga, 89.20.Ff, 89.40.-a 



I. INTRODUCTION 



In daily socio-economic network systems, many com- 
modities, passengers, or information fragments (ab- 
stractly referred to as packets in this paper) are deliv- 
ered from one place to another place. In general, it is 
expected to send and receive packets as quickly as possi- 
ble without encountering traffic congestion which would 
force packets to wait at nodes. Even if packets are con- 
centrated (or condensed in term of statistical physics) on 
just a few nodes in some parts of a network, this situation 
may cause congestion over the network. Thus, one of the 
important issues is a routing strategy: how to select a for- 
warding node in the neighbors of the resident node of a 
packet. Since the efficiency of transportation or commu- 
nication depends on not only routing strategy, but also 
on network topology, we should consider a realistic prob- 
lem setting for the routing and the topology. In addition, 
considering the interaction of accumulated packets in a 
buffer (called queue) at each node is necessary, since it 
crucially affects the traffic flow on a network. This paper 
discusses traffic properties that depend on routing strate- 
gies in the interaction of packets on a realistic network 
topology. 

During this decade, a new research field, complex net- 
work science, has been created [1], since a common topo- 
logical structure called scale-free (SF) was found to exist 
in many real networks such as the Internet, the World- 
Wide- Web, power grids, airline networks, social collab- 
oration networks, human sexual contacts, and biological 
metabolic pathways, etc. The SF structure is character- 
ized by the property that the distribution of degree (the 
number of connections to a node) follows a power-law. 
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In other words, the network consists of many nodes with 
low degrees and a few hubs with high degrees. Moreover, 
a SF network naturally emerges in social acquaintance 
relationships or peer-to-peer communications, and it has 
short paths compared with large network size (the total 
number of nodes). 

In a realistic SF network, a routing path is shortened 
by passing through hubs for the forwarding of a packet. 
However, many packets may be concentrated at the hubs 
[2], in whose queues the packets are cumulatively stored 
if they arrive in a quantity larger than the processing 
limit for forwarding. In such a situation, there exists a 
trade-off between delivering packets on a short path and 
avoiding traffic congestion. To improve communication 
efficiency even in a dynamic environment for a wireless or 
ad hoc network, various routing schemes are being devel- 
oped. Because, in an ad hoc network, many nodes (such 
as base stations or communication sites) and connections 
between them are likely to change over time, then global 
information, e.g., a routing table in the Internet, cannot 
be applied. In early work, some decentralized routing de- 
veloped to reduce energy consumption in sensor or mobile 
networks, however they lead to the failure of guaranteed 
delivery [.':!]; in the flooding algorithm, multiple redun- 
dant copies of a message are sent and cause congestion, 
while greedy and compass routing may occasionally fall 
into inflnite loops or into a dead end. Thus, we focus 
on stochastic routing methods using only local informa- 
tion with respect to the resident node of a packet and 
to the connected neighbors on a path, due to their sim- 
plicity and power. If the terminal node of a packet is 
included in the neighbors of the resident node, then the 
packet is deterministically sent to the terminal in order 
to guarantee reachability unless the connectivity is bro- 
ken. We call this neighbor search (or n-scarch) only at 
the last step to the terminal. Note that the development 
of stochastic routing methods is still in progress depend- 
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ing on device and information processing technologies. 
Although we suppose a mixture of wireless and wired 
communications in future networks, first of all, we aim 
to understand the relations for traffic properties between 
fixed network topology and routing methods. Here, it 
is usually assumed for simplicity that each node has the 
same performance in the forwarding of packets. 

The enhancement of performance in forwarding is also 
reasonable [4, •'>]. For example, in the Internet or airline 
networks, an important facility with many connections 
has high performance; more packets or flights are pro- 
cessed as the incoming flux of communication requests 
or of passengers increases. Such higher performance at 
more busy nodes is effectively applied in a zero-range 
process (ZRP) [6, 7], which tends to distribute packets 
throughout the whole network, instead of avoiding paths 
through hubs. 

In this paper, from the random walk version, we extend 
the ZRP on a SF network to the degree-dependent hop- 
ping rule, which parametrically controls routing strate- 
gies in the trade-off between the selection of a short path 
passing through hubs, and the avoidance of hubs at which 
many packets are condensed and wait in queues for a long 
time. The above models in statistical physics are useful 
for the analysis of complex phenomena involved in tran- 
sition form a free-flow phase to a congestion phase, or 
the opposite transition. Most research on traffic conges- 
tion rely on numerical simulations. Although a few other 
theoretical analyses based on the mean-field approxima- 
tion [s] have been done, we can not compare them to the 
ZRP, simply because of different problem settings. Thus, 
we consider a combination of settings as modified traffic 
models, and numerically investigate them. 

The organization of this paper is as follows. In Sec. II, 
we briefly review the related models in recent studies. In 
Sec. Ill, we introduce a stochastic packet transfer model 
defined by the ZRP, in which all packets persist without 
any generation and removals. Then, in terms of funda- 
mental properties, we approximately analyze the station- 
ary probability of incoming packets at each node on SF 
networks, and derive the phase transition for condensa- 
tion of packets at hubs in the degree-dependent hopping 
rule. In Sec. IV, we numerically confirm the phase tran- 
sition, and as a new result, show the trade-off between a 
detour wandering path and long wait at hubs. Moreover, 
we discuss the congestion phenomenon in the modified 
traffic models with packet generation. In Sec. V, we 
summarize these results and describe some issues for fur- 
ther research. 



II. RELATED WORK 

Many traffic models have been proposed in various 
problem settings for routing, node capacity, and packet 
generation. They are summarized in Table I. The basic 
processes for packet transfer consist of the selection of 
a forwarding (coming-in) node and the jumping-out of 



packets in the node capacity. In these models (including 
our model discussed later), a significant issue commonly 
arises from the trade-off between delivering packets on 
a shorter path and avoiding the congestion caused by a 
concentration of packets on a few nodes such as hubs. 

In the typical models, a forwarding node k is cho- 
sen with probability either K^/^-^j^,Kj [9-11] or 

(™(fc) + l)~''/I]jgAA.("*0) + l)"'' [12]. Here, a and /3 are 
real parameters, Kk and TO(fc) denote the degree and the 
dynamically-occupied queue length by packets at node k 
in the connected neighbors Afi to the resident node i of 
the packet. These methods are not based on a random 
walk (selecting a forwarding node uniformly at random 
among the neighbors), but on the extensions (including 
the uniformly random one at a, /3 = 0) called preferen- 
tial and congestion-aware walks, respectively. Note that 
a > leads to a short path passing through hubs, and 
that a < and /3 > lead to the avoidance of hubs and 
congested nodes with large rn(^k)- In stochastic routing 
methods, instead of using the shortest path, the optimal 
values a = —I and /3 = 1 for maximizing the generation 
rate of packets in a free-flow regime have been obtained 
by numerical simulations [10, 11]. A correlation between 
the congestion at node level and a betweenness centrality 
measure was suggested [12]. 

Other routing schemes [8, 13, 14] have also been con- 
sidered, taking into account lengths of both the routing 
path and of the queue. In a deterministic model [14], 
a forwarding node k is chosen among neighbors Afi by 
minimizing the quantity hdk + (1 — h)m(^k^ with a weight 
< h < 1, dk denoting the distance from k to the ter- 
minal node. Since we must solve the optimization prob- 
lems, these models [11, 14] are not suitable for wireless or 
ad hoc communication networks. Thus, stochastic rout- 
ing methods using only local information are potentially 
promising. In a stochastic model [8], £ A/i is chosen 
at random, and a packet at the top of its queue is sent 
with probability 1 — ri(m(^k)) or refused with probability 
r^im^k)) a-s a nondecrcasing function of the queue length 
m(fe). This model is simplified by the assumption of a 
constant arrival rate of packets, for analyzing the critical 
point of traffic congestion in a mean- field equation [8] . 

With a different processing power at each node [10], it 
has also been considered that the node capacity Cj is pro- 
portional to its degree Ki, therefore more packets jump 
out from a node as the degree becomes larger. On the 
other hand, in the ZRP [6, 7, 15], the forwarding capac- 
ity at a node depends on the number of m^i) defined as 
a queue length occupied by packets at node i. The ZRP 
is a solvable theoretical model for traffic dynamics. In 
particular, in the ZRP with a random walk at a = 0, 
the phase transition between condensation of packets at 
hubs and uncondensation on SF networks has been de- 
rived [(), 7]. For a > 0, a similar phase transition has 
been analyzed in the mean-field approximation [Hi]. 

In the next two sections, based on a straightforward 
approach introduced in Refs. [G, 7], we derive the phase 
transition in the ZRP on SF networks with the degree- 
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dependent hopping rule for both a > and a < 0, in- 
spired by preferential [9-11] and congestion-aware walks. 
Although the rule is not identical to the congestion-aware 
routing scheme [12] based on occupied queue length TO(fc), 
a < corresponds to avoiding hubs with large degrees, 
where many packets tend to be concentrated. Further- 
more, we study the traffic properties in the case with 
neighbor search into a terminal node at the last step. 

TABLE I: Recent traffic models for complex networks. When 
a packet is forwarded from a current node i to k £ Mi, the 
neighboring node k is chosen with probability life or by min- 
imizing an objective function on a routing path. The node 
capacity c; is defined by the number of simultaneously trans- 
ferable packets from each node i. Here, K{xi) denotes the 
degree of node xi, 0{x) is the step function, and m* is a 
threshold, /3 > 0, and < ?? < 1. 



In this routing rule related to the ZRP, the total num- 
ber M of packets is constant at any time, and each node 
performs a stochastic local search as follows: at each time 
step, a packet jumps out of a node i stochastically at a 
given rate qi(m(^^-^) as a function of m^j-), and then comes 
into one of the neighboring nodes k g A/i chosen with 
probability / Y^k' ej\f-^k' ■ When a packet is trans- 
ferred from i to fc, the queue length mj-^j) is increased 
by one unit, and mj-j) is simultaneously decreased. The 
above processes are conceptual, the relaxation dynam- 
ics for the simulation of packet transfer is necessary, and 
described in Sec. IV. The effect of the deterministic 
neighbor search into a terminal node at the last step will 
also be discussed later. 



B. Stationary probability in Q-random walks 
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III. PACKET TRANSFER MODEL 
A. Routing rule 

Consider a system of M interacting packets on a net- 
work of N nodes with a density p = M/N, M = 
SiLi '^(i)i where m^j) > denotes the occupation num- 
ber of packets in the queue at each node 1 < i < N . For 
simplicity, we assume that the queue length is not lim- 
ited, and that the order of stored packets is ignored. If 
there is a limitation on the queue length, it may become 
necessary to discuss a cascading problem [17-19] whose 
dynamics are very complicated. 

To investigate the predicted properties from the the- 
oretical analysis in the ZRP [G, 7, 15], we use the same 
problem setting, such that all packets persist on paths 
in the network without generation or removals of pack- 
ets. For more complex situations with packet generation, 
some models modified by adding different routing rules 
will be investigated in subsection 4.2. 



Before investigating the influence of dynamic queue 
lengths on traffic properties, we consider the stationary 
probability of incoming packets at each node. As men- 
tioned in the previous subsection, we extend a random 
walk routing [15, 2i)], in which a walker (packet) chooses 
a node uniformly at random among the neighbors of the 
current node on a path, to a degree-dependent routing. 
We call it a-random walk. 

As in Ref. [20], for the probability Pij of finding the 
walker at node j € Afk from node i through the inter- 
mediate nodes k G Mi at time t, the master equation 
is 



P^,{t + 1) = 



E 

k 



E 



P.k{t). 



(1) 



By iterating Eq. (1), an explicit expression for the transi- 
tion probability Pij to go from node i through ji, . . . , jt-i 
to i in t steps follows as 
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(2) 



where the sum is taken over the connected 

paths between nodes i and j, as shown in Fig. 1. In 
the opposite directions of the same paths, the transition 
probability Pji follows as 
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We assume the network to be uncorrelated: there is no 
correlation between two degrees of any connected nodes, 
so that K'^, = K, ■ (if") and E,' = ■ (if") 

hold by using the mean value {K") of the a-power of 
degree. By comparing the expression of Pij in Eq. (2) 
with that of Pji in Eq. (3), we obtain the equivalent 
relation 
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Thus, for any source i and step t, Pij is proportional to 
ifj^", while Pji is proportional to K]'^" . Consequently, 
the stationary solution P°° of the incoming probability 

Pij is proportional to K^^"' . This form, related only to 
the degree of the forwarding node j, is suitable for the 
theoretical analysis of the ZRP presented in the next sub- 
section and in the Appendix. The approximative solution 
is also obtained from a different approach to the master 
equation [ID] under the same assumption of uncorrelated 
networks. 




FIG. 1: Combination of consecutive nodes on the route from 
node i to j. The dashed rectangles correspond to the sum of 
Eqs. (l)-(3). 



Figure 2 shows the stationary solution P°° cx Kj and 
the exponent /? « 1 + a > for the SF networks (gen- 
erated by the BA: Barabasi- Albert model [1]). Table II 
also shows that the estimated exponents /3 are consistent 
in the cases both with and without neighbor search into 
a terminal node at the last step. Here, the terminal is 
chosen from all nodes uniformly at random. After ar- 
riving at the terminal, the packet is restarted (resent) 
from the node to a new terminal, in order to maintain 
the persistency of packets in the ZRP. 



C. Phase transition in the ZRP 

We discuss the phase transition between condensation 
and uncondensation in the ZRP with degree-dependent 
hopping rule of packets. For the configuration 

of occupation at each node, the stationary solution of a 
single packet is given by the factorized form 

P(m(i),...,m(,;),...,m(Ar)) = ^HfLiMm^^)), (4) 




K 

FIG. 2: Probability of incoming at a node with degree K in 
the Q-random walk of M — 1000 independent packets (equiva- 
lent to the case of 5 = 1) through 1000 rounds. The lines guide 
the estimated slopes in the left column (without n-search) in 
Table II, and the plotted marks show the probability obtained 
in a BA network with iV = 1000 and (K) = 10. 



TABLE II: Numerically estimated exponent /? ~ 1-fa for the 
lines from the plotted marks in Fig. 2 by using the mean- 
square-error method for the independent walks of M = 1000 
packets. These values give the averaged slopes for Pj° oc Kj 
in 100 BA networks with iV = 1000 and (K) = 10. 
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2.143 
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0.0 
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1.027 


-0.5 
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-1.0 
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0.019 



where Z is a normalization factor and 

for an integer mj^) > and /i(0) — 1. We consider a 
function qi{uj) ~ ui^ with a parameter < 5 < 1 for 
forwarding performance. This means that a node works 
harder for transfer, as it has more packets in the queue 
with larger uj and 5. 

Using the probability distribution in Eq. (4) and 
the stationary probability P^ oc , we can calculate 
the mean value {m^i^) at each node. Here, {m(i)) '= 
'n^=o^Pi^'^{i) ~ ^) defined by using the probabil- 
ity distribution Pi{'m(^i)) = uj) of the number of packets 
occupying node i. 

Ptim{r) =^) = /»(^)nj¥»/i("^(j))' (6) 
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where 



the sum J2* 



E 



is taken in combination 
TO(j;_i_i), . . . , m(7v)} on the constraint 
= M — uj as shown in Fig. 3, and Zi is 



a normahzation factor. It is difficult to directly solve 
the normalization factor Zi or Z. Thus, introducing a 
fugacity variable z [21], the mean value is given by a 
generating function 

E^^^"/.H dlnF.jz) 



dz 



where the second term in the right-hand side of Eq. (7) 
is due to the definition 



m(|| = omega: 

a fixed length 

of occupied queue 



m(,)+... + m(i_,| + m,i^,) + ...+ m,fj)= M - omega 

FIG. 3: Illustration of the queue length occupied by packets 
at each node. For the sum in Eq. (6), there are many 
combinations of {mji), . . . , m(i_i), m(i_|_i), . . . , m(;v)} satisfy- 
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From the definition (5), <7i(aj) — oj^ , and P°° oc iff, we 
have 



(uj \Y ' 

-0 ^ ^ 



(8) 



=1 (^) 



in Eq. (5). The 



because 

fugacity z should be determined from the self-consistency 
equation p = (""^(i))/-^ ^ function of z. Note 

that M ~ is i'h*^ total number of packets in 

a network of N nodes, and that the density p is constant 
at A^, M ^ oo. 

In this paper, we consider a SF network whose degree 
distribution follows a power-law P(K) ^ . Accord- 
ing to the performances of jumping-out, we classify the 

following cases (A) 5 > 5^ =^ /3/(7 - 1), (B) 8 = 5^, (C) 
(5 < (5c, and (D) (5 = for a critical value 5c of the phase 
transition from uncondensation to condensation of pack- 
ets, or the opposite transition. Since the derivation is the 
same as in Refs. [G, 7] at a = 0, except with a slight mod- 
ification for a general value of a (and the corresponding 
/3 K, 1 + a), we briefly review it in the Appendix. We 
summarize the generalized results for < |q;| < 1 in Ta- 
ble III (from that for a = 0). In the cases (B) and (C), 
many packets condensate at nodes with degree K > K^, 
because of the larger exponent /3/(5 > /3. Note that the 
queue length occupied by packets is rewritten as . 
taking into account the dependence on the degree Ki. 



IV. SIMULATION 

In subsection 4.1, we numerically investigate the basic 
properties of packet transfer in the ZRP. In subsection 
4.2, we further discuss congestion phenomena with packet 
generation. 



TABLE III: Scaling of the crossover degree Tfc, the mean 
occupation number m/f at a node with degree K and muub 
at the hub with the maximum degree Kmax for the cases: (A) 
S > 5c, (B) 5 = Sc, (C) < S < Sc, and (D) S = 0. A blank 
denotes no correspondence. 
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Traffic properties for a-random walks in the 
ZRP 



We have performed simulations for M = 1000 packets 
on SF networks generated by the BA model [1] with a size 

= 1000, an average degree (A') = 10, and an exponent 
7 = 3 of P{K) ~ K~'^ . From the initial state in which 
one packet is set on each node, the following processes 
are repeated as the relaxation of the ZRP [7] from node- 
based dynamics to particle-based dynamics [21]. At each 
time step, a packet is selected at random. With prob- 
ability (7i(TO(j))/m(i) = "^(i)^' packet jumps out of 
its resident node i, and hops to one of the neighboring 
nodes j G N^ with probability K'^IY.yaMi Other- 
wise, the selected packet does not move with probability 
\ — qi {jn(i-) ) . The time is measured as a unit of one 
round (Monte Carlo sweep) consisting of M trials of the 
random selection of a packet. 

Figure 4 shows the mean occupation number {ttik) 
of packets at a node with degree K in the cases both 
with/ without neighbor search (n-search), while the top 
of Fig. 4 shows the predicted piecewise linear behav- 
ior [IG] for the crossover degree shown in Table IV; 
the steeper line indicates condensation at the nodes with 
high degrees, and the bottom of Fig. 4 shows that con- 
densation is suppressed by a more gentle line. As shown 
in the insets, the marks slightly deviate from a line, es- 
pecially in the head and the tail, because of the effect of 
n-search into a terminal node at the last step. In these 
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TABLE IV: Classification of tlie cases in Table III for the 
combinations of a and 5. The values in the mid-columns are 
Kc = Kl^aJ^" , for 5c = 13/2 in the case (C), corresponding to 
the maximum, average, and minimum of Kmax ~ 189, 132, 98 
in each vertical triplet from the top to the bottom. These 
results are obtained in 100 realizations of BA networks with 
N = 1000 and {K) = 10. 
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cases at the top and the bottom of Fig. 4, the critical 
values of <5c « (1 + a)/{j - 1) are 0.75 > 6 and 0.25 < S, 
respectively. Thus, condensation of packets at hubs can 
be avoided as the performance of transfer is enhanced at 
a large S, although the transition depends on the prob- 
ability of incoming packets according to the value of a; 
a negative a induces a nearly homogeneous visiting of 
nodes, while a positive a induces a heterogeneously bi- 
ased visiting of the nodes vifith high degrees. We note that 
the uncondensed phase is maintained in a wide range of 
5 > Sc ^ {1 + ct)/{j — 1) for a < 0, as shown in the case 
(A) of Table IV. 

Next, in order to study the traffic properties, such as 
the travel time on the routing path, we consider packet 
dynamics in a realistic situation with n-scarch into a ter- 
minal node at the last step. If the terminal is included 
among the neighbors of the resident node for a randomly 
selected packet, then it is deterministically forwarded 
to the terminal node, taking into account the reach- 
able chance, otherwise it is stochastically forwarded to 
a neighboring node j with probability K" / Kp. 
The n-search is practically effective and necessary in or- 
der to reach a terminal node. Remember that the esti- 
mated values of /3 are similar to both with/without n- 
search, as shown in Table II. Therefore, the behavior 
of the mean occupation number (tok) is similar to that 
shown in Fig. 4 and the inset. In the following discus- 
sions, if we leave out n-search, then the optimal parame- 
ter values of a and S may be changed for low-latency de- 
livery, although the difference is probably small from the 
above similarity. By selecting hubs for a > 0, a packet 
is terminated with higher probability even in only the 
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Case(C): a = 0.5, S = 0.2 
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Case(A): a = -0.5, S = 0.8 

FIG. 4: Typical results for the mean occupation number {uik) 
of packets at a node with degree K in the routing without n- 
search. Inset: the cases with n-search. The circle and cross 
marks correspond to the results in 1000 and 100 rounds, re- 
spectively. The dashed lines show the slopes /3 and [3/5 in 
Table III. Note that condensation of packets at the nodes 
with high degrees occurs in Case(C). These results are ob- 
tained from the averages of 100 samples of packet transfer on 
a BA network whose maximum degree is the closest to the 
average value Kmax ~ 132 in the 100 realizations. 



stochastic forwarding, because a terminal node is highly 
likely to be connected to some hubs. At the same time, 
this leads to congestion at hubs. However, in the case 
without n-search, a packet may wander for a very long 
time, which is not bounded a priori for the simulation of 
the packet transfer. It is intractable due to huge compu- 
tations. Thus, we focus on the case with n-search. 

We investigate the traffic properties for reachability of 
packets, the number of hops, the travel time Tq including 
the wait time of a trapped packet in queues, and the sum 
Tw of wait times on a path until arrival at the terminal. 
Note that the travel time is Ta = T^, + (time to one 
hop) X (the number of hops). We also define the averaged 
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FIG. 5: Convergence properties in the high-performance 
regime at 5 = 0.8 in the observed interval [rounds] after the 
transient state of (mif). Inset: results in the low-performance 
regime at 5 = 0.0. These results are obtained from the av- 
erages of 100 samples of packet transfers on a BA network 
whose maximum degree is the closest to the average value 
Kmax = 132 in the 100 realizations. 



wait time T^/N^ per node, where is the number 
of trappings in queues at the nodes on a routing path. 
These measures are cumulatively counted in the observed 
interval after 1000 rounds, and averaged over 100 samples 
of this simulation. Here, we discarded the initial 1000 
rounds as a transient before the stationary state of TO(i) 
is reached. 

Figure 5 shows, from the top to the bottom, the con- 
vergence of reachability, the mean number of restarted 
packets per round, and the mean travel time in the high- 
performance regime aX 5 — 0.8. The inset shows a 
slightly slow convergence in the low-performance regime 




FIG. 6: The mean travel time (top) and the mean number of 
hops (bottom). The A, -f , *, x, and V marks correspond to 
a = 1.0, 0.5, 0.0, —0.5, and —1.0, respectively. The simulation 
conditions are the same as in Fig. 5. 



at 5 = 0.0. For other measures, similar convergence prop- 
erties are obtained. We compare these curves in terms of 
the values of a; they shift up from a = —1.0 to a = 0.5 in 
the case (A) 5 > 5c, while down from a = 0.5 to a = 1.0 
in the case (C) 5 < Sc- This non-monotonic dependence 
on a is related to the condensation transition, since the 
change between 6 > 6c and 6 < Sc occurs at the critical 
value a = 0.6 for the equivalence 6c ~ {1 + a)/{'-f — 1) 
(5 = 0.8. Thus, it appears when a is greater or less 
than 0.5 in Fig. 5. In the following, we briefly explain 
each property: reachability, the number of restarts, and 
{To). The reachability of the restarted packets is around 
0.99 on similar curves for all values of a at ^ = 0.8, 
while these curves separate in increasing order of a at 
^ = 0.0 (see the inset at the top of Fig. 5). Note that 
the number of restarted packets cumulatively increases 
as the observed interval is longer, although the rate per 
round is constantly around 3 ~ 5 in the ordering from 
a = 0, a = ±0.5, to a = ±1 for S = 0.8, and less than 
1 in the ordering from a = —1.0 to a = 1.0 for 6 ~ 0.0, 
as shown in the middle of Fig. 5 and in the inset. The 
maximum and the minimum lines for the mean travel 
time (Tq) are obtained at a = — 1 and a = 0, respec- 
tively, for 6 = 0.8. However, all of the curves shift up 
for S = 0.0 (see the inset at the bottom of Fig. 5), and 
(To) is longer as the value of a increases, because of the 
waiting at high- degree nodes. 

We further investigate the above traffic properties, es- 
pecially for the forwarding of packets with more detailed 
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values of 6 for 30000 rounds in the quasi-convergent state. 
As shown in Fig. 6, the mean travel time (To) decreases 
as the value of S increases with higher forwarding per- 
formance, because the wait time trapped in a queue de- 
creases on average. In particular, packets tend to be 
trapped at hubs for a long time when a > 0, and then 
(Ta) is longer. In contrast, the number of hops increases 
on average as the value of i5 increases, because some 
longer paths are included in higher reachability. The path 
length counted by hops tends to be short through hubs 
when a > 0, however it tends to be long on a wandering 
path when a < 0. Therefore the number of hops increases 
in decreasing order of a. Note that the number of hops 
is very small compared to (Ta), which is dominated by 
the mean wait time (T^,) ss (Ta) on a path. 

As shown in Fig. 7, the mean number (Nw) of trap- 
pings at nodes on a path increases when a > and de- 
creases when a < as the value of S increases with higher 
forwarding performance. This up-down phenomenon 
may be caused by a trade-off between the avoidance 
of trapping through sufhciently high forwarding perfor- 
mance at a node and the inclusion of longer paths with 
high reachability. The mean wait time (Tw/Nw) per node 
decreases as the value of 6 increases, and is longer in in- 
creasing order of a because of the long wait time at hubs. 

In summary, the wandering path for a < better re- 
duces the mean travel time of a packet with high reacha- 
bility in the low-performance regime at a small 6, while in 
the high-performance regime at a large S, the difference 
between a > and a < is small, neither the wandering 
long path with short wait trapped at nodes (a = — 1), 
nor the short hopping path with long wait trapped at 
hubs (a ~ 1) is advisable. Thus, a uniformly random 
walk (a = 0) yields slightly better performance. 

B. Congestion phenomenon 

This subsection discusses the congestion phenomenon 
when packets are randomly generated at each node at a 
constant rate p, and removed at the terminal nodes (not 
restarted within the persistency) . In order to reduce the 
computational load, the packet dynamics starts from the 
initial state: there arc no other packets than the ones 
that are generated. We compare the phenomenon in our 
traffic model based on the ZRP with that in the following 
modifications related to the traffic- aware routing [1 4] at 
a = [n] . Table V summarizes a combination of the basic 
processes: with or without (Yes or No) the refusal of for- 
warding, n-search, and a constant arrival with probability 
fi. The other processes for choosing a forwarding node 
j with probability oc K°' and for jumping-out a packet 
from its resident node i at the rate TO*j-| are common. 

Mod 1: With probability r]{m(^j)) ~ fj<d{ra(^j-^ — m*), the 
selected packet is not transfered to j G Afi, but is 
refused at node i, where Q{x) is the step function, 
m* is a threshold, and a parameter < ^ < 1. 
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FIG. 7: The mean number of trappings at nodes (top) and 
the mean wait time per node (bottom). The A, +, *, x, 
and V marks correspond to a = 1.0,0.5,0.0, —0.5, and —1.0, 
respectively. The simulation conditions are the same as in 
Fig. 5. 

Here, node j is chosen deterministically by n-search 
if j is the terminal node, otherwise it is chosen with 
a probability proportional to K". 

Mod 2: Moreover, instead of n-search, the selected 
packet is removed with probability /i, or it enters 
the queue with probability 1 — /i. 

Mod 3: In Mod 2 , there is no refusal process (^ = 0). 

In Mod 1, the randomly-selected packet from the queue 
does not leave node i with a constant probability fj if 
the occupation number Tj^j) of packets is greater than a 
threshold to*, while in Mod 2, after passing this refusal 
check, it is removed with a constant probability ^. Note 
that constant arrival was assumed to theoretically predict 
the critical point of traffic congestion in the mean-field 
approximation as iV H> oo [8]. With packet generation, 
we can perform the ZRP as an extension of the model in 
Rcf. [s]. In particular, the case oi S — corresponds to 
node capacity Ci ~ 1 for all nodes i; only one packet is 
transferable from a node at each time. 

For a variable generation rate p, we investigate the 
appearance of congestion by using the order parameter 
[8, 14] 

M{t + T)-M{t) 
op = hm , 

t^oo Tply 
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TABLE V: Modified traffic models. 



Refusal 


n-search | const . arrival 


No 77 = 


ZRP 


Mod 3 


Yes r? > 


Mod 1 


Mod 2 



where M{t) denotes the sum of existing packets in queues 
over network (praetically for a large t), and r is the ob- 
served interval. The value of op represents level of eon- 
gestion, e.g. op « indicates a free-flow regime, while 
op « 1 indicates a congested regime. 

In the following, we set fj = 0.7 and m* = 5 for the re- 
fusal process. As shown in Fig. 8, in the cases in which n- 
search takes place, the value of op rapidly grows with the 
increasing of the generation rate p, since the removal of 
a packet arriving at the terminal node is rare, especially 
for a small 6. Here, the marks and color lines indicates 
different values of S and a. By comparing the correspond- 
ing curves of the same color and marks in the top and 
the bottom figures, we notice that the ZRP suppresses 
congestion in smaller ops than Mod 1, in particular, for 
6 = 0.5,0.8 (blue and magenta lines). At the top of Fig. 
8, the magenta lines indicate the existence of a free-flow 
regime around a small p for S = 0.8 as high forwarding 
performance. Thus, the refusal process does not work 
effectively in the cases with n-search, at least for this pa- 
rameter set. At both the top and the bottom of Fig. 8, 
the difference for the same color lines with three marks 
corresponding to a = ±1,0 is small, except for S = 0.8 
in Mod 1 (magenta lines at the bottom). However, the 
curves shift down as i5 becomes larger; in other words, 
the congestion is suppressed by higher forwarding per- 
formance. Note that a uniformly random walk at a = 
(asterisk marks) yields better performance for each value 
of (5. 

When a constant arrival with probability ^ is applied 
instead of the realistic n-search, the behavior changes. 
Figure 9 shows that the refusal process works effectively, 
since Mod 2 with the refusal process (at the bottom) has 
smaller op than Mod 3 without the refusal process (at the 
top). The gap between three marks for the lines of each 
color for Mod 2 at the bottom of Fig. 9 resembles to that 
at the top of Fig. 8, however the gap appears remarkably 
at (5 = 0.8 (magenta lines) for Mod 3 in the top of Fig. 9. 
Although n-search leads to a low arrival rate (lower than 
/i = 0.01) as shown in Fig. 5, and reachability is not 100% 
even in the case of persistent packets after restarting, in 
the meaning of smaller op, the ZRP is better than other 
models, by comparison with the corresponding curves in 
Figs. 8 and 9. 

We consider the other stochastic routing method in 
which a forwarding node fc e A/i is chosen with probabil- 
ity 

Hfc cx Kkim^k) + l)-^ (9) 
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FIG. 8: The values of op for the ZRP (top) and Mod 1 (bot- 
tom) using n-search. The open A, *, and v marks correspond 
to Q = 1.0,0.0, and —1.0, respectively. The red, green, blue, 
and magenta lines correspond to 5 = 0.0, 0.2, 0.5, and 0.8, 
respectively. 

where 13 — 3 yields the maximum generation rate in a 
free-flow regime [13]. As shown in Fig. 10, in this optimal 
case, the behavior is similar to that in Mod 3 without the 
refusal process at the top of Fig. 9, although it has bet- 
ter performance than Mod 3. There is a free- flow regime 
in the case when 6 — 0.8 (the magenta line with filled 
upward-pointing triangle marks) at a constant arrival. In 
this method, a forwarding node is dynamically selected in 
a balance between reducing distance by passing through 
large degree nodes, and avoiding congestion. Thus, a 
further improvement may be potentially expected in the 
tuning of the balance for a-random walks and other mod- 
ifications. 



V. CONCLUSION 

For a SF network, whose topology is found in many 
real systems, we have studied extensions of the ZRP [G, 
7, 15, 20] which controls both the routing strategies in the 
preferential [9-11] and congestion-aware [12] walks, and 
the node performance for packet transfers. Under the 
assumption of persistent packets, we have approximately 
analyzed the phase transition between condensation of 
packets at hubs and uncondensation on SF networks by 
another straightforward approach [(i, 7] instead of the 
mean-field approximation in the preferential walk for a > 
[16]. In particular, we have found that uncondensation 
is maintained in a wide rage of J > i5c ~ (1 + Q^)/(7 ~ 1) 
for a < 0. 
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TABLE VI: Qualitative summary of the traffic properties. 
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FIG. 9: The values of op for Mod 3 (top) and Mod 2 (bottom) 
at a constant arrival with probability fj, = 0.01, instead of 
using n-search. The filled A, x, and y marks correspond to 
a = 1.0,0.0, and —1.0, respectively. The red, green, blue, 
and magenta lines correspond to 5 = 0.0, 0.2, 0.5, and 0.8, 
respectively. 
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FIG. 10: The values of op for a trafBc-aware routing [l.i] in 
which a forwarding node k is chosen by applying Eq.(9) in 
the cases of n-search (open marks), and of a constant arrival 
with probability — 0.01 (filled marks). The red, green, 
blue, and magenta lines correspond to 5 = 0.0, 0.2, 0.5, and 
0.8, respectively. 



Moreover, we have numerically investigated the traf- 
fic properties when the practical n-search into a termi- 
nal node at the last step takes place. The phase tran- 
sition has been consistently observed in the cases both 
with/without n-search. The obtained results are summa- 
rized in Table VL We conclude that the wandering path 
for a < better reduces mean travel time of a packet 
with high reachability in the low-performance regime at a 
small 6, while in the high-performance regime at a large S, 
neither the wandering long path with short wait trapped 



Mean Value 


a > 


Q < 


for larger 5 


Reachability 


low 


high 


increased 


Num. of Hops 


small 


large 


increased 


Travel Time (T^) 


long 


short 


decreased \j 


One Wait {T^-)/{N^) 


long 


short 


decreased \, 


Num. of Wait (TV,^) 


small 




increased 






large 


decreased \j 


Characteristics 


condensation 


uniformly 




of a Path 


at hubs 


wandering 





at nodes {a = —1), nor the short hopping path with long 
wait trapped at hubs (a = 1) is advisable. A uniformly 
random walk (a = 0) yields slightly better performance. 
This optimality at a = in a high-performance regime is 
consistent with the results obtained for the critical gen- 
eration rate in other traffic model [ID] on a SF network 
with node capacity (corresponding to the forwarding per- 
formance) proportional to its degree. However, in the 
details for our traffic model, the optimal case depends 
on a combination of the values of a and d related to the 
condensation transition at Sc- We emphasize that, for 
high reachability, small number of hops, and short travel 
or wait time, such traffic properties summarized in Table 
VI as the trade-off between a detour wandering path and 
long wait at hubs cannot be obtained from only the theo- 
retically predicted phase transition between condensation 
and uncondensation in Ref. [6, 7, 1-5, 16]. Concerning 
the fundamental traffic properties, we have investigated 
the congestion phenomenon with packet generation, com- 
pared with other models related to the traffic-aware rout- 
ing [(S, 1:], 14]. We suggest that the high-forwarding per- 
formance at a large 6 is more dominant than the refusal 
process in order to suppress congestion in a small op, and 
that the difference of op is small when the values of a are 
varied (the case of a = is slightly better for n-search). 
The above results only show qualitative tendencies. In 
further research, more complex and quantitative proper- 
ties should be carefully discussed for many combinations 
of parameters a, i5, /i, fj, to*, etc., although such simula- 
tions may be intractable due to huge computation load 
and memory consumption for the processing of millions 
of packets in very long iterations. 

This direction of research on stochastic routing may be 
the first step to reveal complex traffic dynamics such as 
the trade-off between the selection of a short path and 
long wait (delay) at some particular nodes related to the 
underlying network structure. We will investigate this 
problem, including the effects of more realistic selections 
of the source and the terminal nodes which depend on 
geographical positions or population density, the queue 
discipline such as FIFO or LIFO, and other topologies 
to develop the optimal routing schemes [22] for advanced 
sensor or ad hoc networks. 
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Appendix 

As a modification to Refs. [G, 7], we briefly review the 
derivation of the results presented in Table III. 

For (5 = 0: case (D), we have a constant jumping rate 
Qiiuj) = 1 independent of the occupation number uj and 



-. Then, from Eq. (7), 



zK': 



is divergent for the hub with the maximum degree Kmax 
as z approaches Zc = l/K^^^^ in the range z < Zc, be- 
cause the denominator becomes zero in the first term 
of the right-hand side, equivalently to the case in which 
node i has the maximum degree in the second term. It is 
convenient to decompose the density into two terms such 
that p = ps + Pn for the hub with the maximum degree 
and for the other nodes. These two terms are given by 



mhub 
N 



Nl 



zK^ 



(A.l) 



-1 zK^ 



K„ 



1 - zKf> 



P{K)dK. (A.2) 



By using the relation ~ X^z^o -f^^- i^-'^)^ 
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derive 



j^PiK)dK 
SzK0Y.t,{zKPyP{K)dK 

1^1=1 \ ^ 01-1+1 J 



lii 



/3 



(=1 



_ [-ln(l-z/f^,J] 
-iV-i xln(l + mft„b)^0, 

for a large size N oo. In the last terms, we apply 



(-7+1) 



AT-i from /^^ P{K)dK - 1/A^, and 



1 °° ; 



In 



1 - zK'„ 



1 + psN = 1 + m/i„6 > 0, 



from Eq. (A.l). Since only ps remains, condensation 
occurs at the hub for any a in case (D) at this uniform 
jumping rate given by (5 = 0. 

In the general case where i5 > 0, the infinite series in 
Eq. (8) does not have a closed form, except for Fi{z) = 
e^^i at (S = 1. We approximate the series as in [(i, 7], 

« -L(24,iff X exp(5(zXf 

The second term in the right-hand side is dominant for 
zK^ > 1. Wc can simply approximate it with a few 
lowest-order terms such as Fi (z) = 1 + z-ftTf 4- 0{{zK^)'^) 
for zK^ < 1. 

From Eq. (7) and d\n Fi{z)/dz for the above two cases 
of Fi{z), the mean occupation number is given by 



rriKi 



zK^ for zK^ < 1 

{zK^f/^ for zK^ > 1. 



Thus, the mean occupation number of a node increases 
monotonically with its degree. Depending on a constant 
magnitude of z, we consider the following two cases (i) 
and (ii). 

(i) We assume that the fugacity z is in such a range 
that zK^ > 1 for all nodes. From rriKi — {'^^iY^^ fo^' 



nodes and the self-consistent equation P = j^'^i "^A'i 
TfZ^/^Y.r{Kt?'\ the solution z = p^/{kf^/^Y is valid 
only when Rf^'^ =^ /^"^°" Kf^'^P{K)dK remains finite 
in the limit of large N. By using P{K) ^ K^'', at TV ^ 
oo, we have 



K^/^P{K)dK 



K„ 



K^"-"'dK^K^£'i-''+'\ 



From the negative exponent, this finite condition imposes 
that 5 > 6c with 



dcf 

Or, = 



/3 



7- 



In this regime, we find . K^^^ for all nodes. The 
occupation number at the hub with the maximum de- 
gree K^nax ^ A^i/(')'-i) scales sublinearly as rrihub ^ 
jV/3/5(7-i) ^ Tv-^c/s ^i^ij ti^e exponent /3/(5(7 - 1) = 
5c/ 5 < 1. Therefore, the extreme condensation of al- 
most all packets at the hub with K„iax is avoided when 
5 > 5c- case (A). 

(ii) We assume that the fugacity z is defined as = 
l/z in the interval Kmin < Kc < K^ax- Then, the 
self-consistent equation becomes p = Pn + Ps = X^i ; 
where 



Pn= J2 = 



K„ 



■f^ I K^P{K)dK, (A.3) 



is the density for the nodes with K < Kc- zK^ < 1, and 



(A.4) 

is the density for the nodes with K > Kc'- zK^ > 1. 
Since the integral part in Eq. (A.3) is smaller than the 
average /3-power of the degree, pn vanishes as K'^ in the 
limit of large TV. In order to have a finite value of ps = p, 
the integral part in Eq. (A.4) should be divergent, which 
yields that 5 < 5c — P / {"^ — 1). This integral should be 
of the same order as Kc which yields 



Kc 



[XnK^axf^'^ for 5 = 5c 
{Krnaxf-'"^ for 5 < 5c. 



in cases (B) and (C), respectively. 



